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Abstract 

We discuss quantum tunneling between classically BPS saturated solitons 
in two-dimensional theories with N = 2 supersymmetry and a compact space 
dimension. Genuine BPS states form shortened multiplets of dimension two. 
In the models we consider there are two degenerate shortened multiplets at 
the classical level, but there is no obstruction to pairing up through quantum 
tunneling. The tunneling amplitude in the imaginary time is described by 
instantons. We find that the instanton is nothing but the 1/4 BPS saturated 
"wall junction," considered previously in the literature in other contexts. Two 
central charges of the superalgebra allow us to calculate the instanton action 
without finding the explicit solution (it is checked, though, numerically, that 
the saturated solution does exist). We present a quantum-mechanical inter- 
pretation of the soliton tunneling. 
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1 Introduction 



Bogomol'nyi-Prasad-Sommerfield (BPS) saturated topologically stable solitons in 
supersymmetric theories are widely discussed at present in connection with the 
brane world scenarios 0,0. In theories with compact spatial dimension(s), two 
distinct degenerate mass solitons which are BPS saturated classically, and to any 
finite order in perturbation theory, can mix nonperturbatively, thus lifting the BPS 
bound. Two shortened super mult iplets pair up with each other combining in one 
full supermultiplet with mass M > \Z\, where Z is the central charge of the super- 
algebra. This phenomenon is an analog (in the soliton sector) of the spontaneous 
breaking of supersymmetry due to instantons in the vacuum sector [[], f| . To the best 
of our knowledge it was first considered in the context of M = 2 two-dimensional 
Wess-Zumino models in 

In this work we address the issue of calculating the shift M — \Z\. In the quasi- 
classical approximation, it is proportional to the tunneling probability which, in 
turn, is determined by instantons. Remarkably, the instanton calculus in this case 
is nothing but an adaptation of the theory of the BPS saturated wall junctions, 
which also received much attention recently || [7], ||, ||, [10], [II], |T2| . In particular, the 
instanton action can be derived from the central charges. The explicit formula for 
the instanton solution is not needed. The only thing we need to know is the very fact 
of its existence. This is in perfect parallel with the standard (non-supersymmetric) 
instantons: once one knows that the self-duality equations have a solution, the 
instanton action is unambiguously fixed in terms of the topological charge. 

In Ref. || it was observed that the soliton mixing, resulting in the loss of the BPS 
saturation, can be described by an effective SUSY quantum mechanics; however, the 
general construction presented there, is not very transparent. Here we reduce the 
construction of Ref. [[J to a simple limiting case which nicely illustrates the essence 
of the phenomenon. 

The organization of the paper is as follows. In Section 2 we formulate the problem 
and elaborate general aspects of the solutions. In Section 3 a specific instructive 
example is considered. Section 4 is devoted to SUSY quantum mechanics. 



2 Formulation of the problem and general results 

Classically BPS saturated soliton supermultiplets which may become degenerate in 
mass with some other supermultiplets and lift the BPS bound because of a non- 
perturbative mixing, is a rather general feature of various theories. Although our 
results are applicable in all cases, we find it convenient to explain the problem in a 
specific setting. 

Consider a two-dimensional M = 2 Wess-Zumino model of one chiral superfield 
$ with the superpotential W($). Any model of this type can be obtained as a two- 
dimensional reduction of the corresponding four- dimensional theory. The geometry 
of the world sheet is a cylinder. As explained in ||, for the existence of the BPS 
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Figure 1: a) The general singularity structure of the quantity dW/d§ in the complex 
$ plane; b) the critical trajectories corresponding to infinite periods. 



solitons it is necessary that W($) is a multi-branch function (otherwise the vanishing 
of the central charge Z cannot be avoided), while dW/d§ must be meromorphic. 
Another necessary (and sufficient) condition for the topologically stable solitons is 
the existence of non contractable cycles in the target space. One of the simplest 
choices is a target space with the topology of a cylinder, possibly with punctured 
points. Then the periods of W are the central charges of the SUSY algebra, 

dW, (1) 

where nc; stands for the i-th non-contractable contour in the target space. 

If there is at least one non-contractable cycle, one can always define $ in such a 
way that dW is periodic 

ciW($ + 2vr) = dW($). (2) 

This particular parametrization is not crucial, and is imposed only for the purpose 
of making the presentation simpler. The poles of dW/d& are assumed to be single 
poles; quadratic and higher order poles can be treated as a limiting case of coinciding 
single poles. A generic singularity structure of dW/d& is depicted in Fig. [T|a), where 
the poles are marked by bold dots. The Kahler potential is taken to be trivial, 

/C($,<5) = 

Consider the cycle Ti, for which 

Z 1 = 2/«('!^V (3, 



If Z\ = \Zi\e lS , then the equation for the static BPS soliton has the form 

dx d§' y ' 
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This equation admits the "integral of motion" , 

/ = Im (e- tS W) , (5) 

i.e. dl/dx = when W and VV are evaluated on the solution of (f|). The existence 
of this integral of motion allows one to find the BPS solution in the general case ||. 
The strategy is as follows. We first ignore that the world sheet is a cylinder with 
period L in the x direction, and solve the BPS equation without posing the condition 
of periodicity, $(x + L) = The solution found in this way is marked by the 

continuous (real) parameter J, and is a periodic function of x with period 

J \ aw J $=$ J ( a; ) 

The period function £(I) is real and positive. It is not difficult to show that there 
exist critical values of / such that £{I) — > oo at / J, i*/ or / f I vn where I*i^ r mark 
the critical trajectories running close to the nearest poles of W from the left and 
from the right (Fig. |I|b). A schematic plot of the function £(I) is presented in Fig. |2|. 

If the circumference of the worldsheet cylinder L > £ m i n , then the equation 
£(I) = L has an even number of solutions. The corresponding values Jj belongs to 
the interval (i*/,i* r ), while $^(x), i = 1,2, . ..,2u are the classical BPS solutions 
satisfying Eq. (f|) and the periodicity condition $^(x + L) = $^(x). 

In the case at hand they have particle interpretation. Altogether, we have 2v 
supermultiplets, each containing two states. All masses are degenerate and equal to 

The BPS nature of the solitons established above at the classical level persists 
to any finite order in perturbation theory (this statement assumes that there is a 
small expansion parameter in the superpotential and/or the Kahler function). Al- 
ternatively, one can say that 1v BPS solitons remain under small deformations of 
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the parameters. This is due to the fact that the number of states in the supermul- 
tiplet is two, while the full Af = 2 supermultiplet contains four states. The BPS 
supermultiplets are shortened. 

It is equally clear, however, that nonperturbatively the BPS saturation of the 
solitons under consideration is lifted, they pair up to form v full supermultiplets 
which lift the BPS bound. This was noted in ||, where arguments were given 
based on the Cecotti-Vafa-Intriligator-Fendly (CVIF) index Our task here is 
to calculate M — \Zi\ in the quasi-classical approximation (which implies of course 
that (M — |Zi|)/|Zi| < 1). 

The BPS saturation is lifted by tunneling. Consider for simplicity the case when 
Eq. (H) has only two solutions, $i t (x) = 4>x{x) and $i 2 (x) = 2 (x)- One can 
construct an interpolating field configuration <f>(t, x) (where t is the Euclidean time) 
such that in the distant past 

<p(t,x) ><Pi(x), (7) 

and in the distant future 

<f>(t,x) t=T/2 ^°° y falx), (8) 

The interpolation is smooth (in particular, <f)(t,x + L) = <f)(t,x) for all t), so that 
the (Euclidean) action A 

A[<P(t,x)]-\Z\T (9) 

is finite. Here \Z\ is the soliton mass in the absence of the tunneling. The quasi- 
classical formalism is applicable provided A — \Z\T ^> 1. One must minimize over all 
interpolating trajectories; once the trajectory (f)o(t,x) corresponding to the minimal 
action is found, one can calculate 

AA min = A[(j> (t,x)}-\Z\T. (10) 

The shift of the soliton masses from the BPS bound is then 

AM = M - \Z\ oc e ~ 2AAmh \ (11) 

where the factor of 2 in the exponent is due to the fermion zero modes. We will 
comment more on this factor in Section 4. 

The central result of the present work is as follows. The determination of the 
minimizing trajectory (f>o(t,x) (in the Euclidean time) reduces to the problem of 
determining the BPS wall junction in the (1 + 2)-dimensional theory, in which the 
(1 + l)-dimensional model under consideration is embedded. The embedding is triv- 
ial. Indeed, if the original model is a (1 + l)-dimensional slice of the four dimensional 
Wess-Zumino model, the one in which we embed is a (1 + 2)-dimensional slice of 
the very same model. From this remark it is clear that the formalism we discuss 
is applicable, generally speaking, in the supersymmetric theories with extended su- 
persymmetry (A/ - = 2 or M = 4). In fact, since the solitons - the "walls" and "wall 
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Figure 3: The energy distribution corresponding to the interpolating configuration; 
the unshaded has the "vacuum" energy density. 



junctions" - are static, we will have to deal with one- and two-dimensional problems, 
respectively. 

The energy distribution for the interpolating configuration <fr (t,x) is schemati- 
cally depicted in Fig. £| It is nothing but an adaptation of the standard four-wall 
junction on the cylinder world sheet. The equation for the BPS wall junction has 
the form 

d( ~ 2 d4> ' 1 ] 

where ( is the complex variable ( = x + it, and = 1/2 (d x — idt). The solution 
to this equation, if it exists, is 1/4 BPS saturated. Equation ( |T2"D was first derived 
in The fact that the solution of Eq. ( |T2"D minimizes the Euclidean action is 



quite obvious. Indeed, 
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+ surface terms. 



The surface terms are unambiguously fixed by the boundary conditions, Eqs. (^) 
and (H), and by the periodicity condition <fi(t, x+L) = <p(t, x). Thus, the construction 
is quite analogous to the instanton self-duality equation in the Yang-Mills theory 
or the two-dimensional 0(3) sigma model, where the surface terms are topological 
charges. 
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In the supersymmetric model under consideration the surface terms are propor- 



tional to two distinct central charges which exist in the superalgebra fLOfl . One 
central charge is related to the junction "spokes". In fact, this was discussed above, 
see Eq. (0). Another central charge is related to the junction "hub". 

In what follows we will assume the circumference of the world sheet cylinder to 
be large. This is sufficient to ensure the applicability of the quasi-classical approxi- 
mation. 

In the quasi-classical limit L — > oo the central charge related to the junction is 
subdominant. From Fig. |^, it is evident that at L — > oo the minimal action is 

AA min = La, (14) 

where a is the tension of the horizontal "wall" , 

<y = 2{h r -h l ). (15) 

The effect of the "hub" is subdominant, it is proportional to L°, and may be of the 
same order as the pre-exponential factors due to the zero modes. Thus, 

M- \Z\ oc e - 4{hr - hl)L . (16) 

In the next section we will consider a concrete model, which seems to present an 
instructive example. In this particular model we calculate for AA min both the linear 
term in L and, for the sake of completeness, the next to leading term associated 
with the (1/2, 1/2) central charge, the "hub". 

The remainder of the paper presents an illustration and elaboration of the above 
general results. 



3 An (instructive) example 



In this section we apply the previous considerations to a specific model which was 
first introduced in ||. We consider a generalized Wess-Zumino model for which 

- 4-7T 

= dW = d$, (17) 

2 — cos $ 

where dW/d<& is a single- valued function derived from the multi- valued superpoten- 
tial 

8n ( r & 

Yv = arctan v 3 tan - 

Vs V 

The model possesses only a run-away vacuum |lm(0)| — > oo. It is stabilized by 
solitons, which at the classical level are solutions to the BPS equation given by 
Eq. (§. 
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The target space has the topology of a cylinder (—00 < lm(0) < +00, — n < 
Re(0) < 7r), with the poles of the scalar potential, 



11,2 



±zlog(2 + V3), (19) 



removed. Each soliton solution belongs to one of three homotopy classes, r 1; T 2 and 
T 3 . Solutions in Ti wind around the target space cylinder, whereas solutions in T 2 
and T 3 wind around the points that are removed from the target space (see Fig. f|). 
Each solution in Ti has a mirror image in the real axis of the complex plane that 
belongs to the same class, except for a real soliton, which is invariant under this 
transformation. Solutions in T 2 are mapped to solutions in T 3 and vice versa. The 
solutions in all homotopy classes have equal period, 

8tt 2 

AW = 7T (20) 

As explained in Section |2|, the constant of motion / = Im (W) (remember, 5 = 
0) may be used to mark all solutions to the BPS equation. The BPS solitons 
in homotopy class T 2 and T 3 are obtained for / £ (—00, J*;) and / £ (I* r ,+oo) 
respectively, with 7^ = —I c and I vr = +I C , and 

8tt , 1 

I c = ^^arctanh^^. (21) 

The classically BPS saturated solitons in homotopy class Ti are obtained for / £ 
(hi, I„). The period function 1(1) was plotted in Ref. |J; this function is positive, it 
is symmetric under reflection in I = 0, and it monotonically increases from £(I) J, 
at / — > — 00 to £ — ► 00 at / t I*i- Between I = I*i and I = I* r , the period function 



7 



reaches a minimum value t — 1 for I — 0, where the classically BPS saturated 
soliton is real. 

For a given circumference L of the compact dimension, the allowed solitons with 
winding number N are obtained from the equation N£(I) = L. The energy of BPS 
saturated solitons is independent of L and given by 

, , 16iV7r 2 , s 

M BPS = \Z\ = (22) 

In this paper we explicitly discuss solitons with winding number N — 1, but the 
results can be trivially extended to arbitrary values of N. In the present model, 
there are two BPS saturated solitons if L < 1, one in homotopy class T 2 and one 
in r 3 . If L > 1 there are four classically BPS saturated solitons, two in homotopy 
class Ti and one each in T 2 and T 3 . 

It turns out that for practical purposes, it is convenient to mark the class Ti 
solitons by the quantity B, the imaginary part of <fi when the real part of is equal 
to 7r. There is a one to one correspondence between B and / that takes the form 



57T 



-arctanh ( —= tanh — , (23) 
so that for the critical solitons B(±I C ) — > ±00, and B — for I — 0. 



1 . B 



3.1 Non-BPS solitons for L < 1. 

Even though there are no BPS solitons in homotopy class Ti when L < 1, this does 
not preclude the existence of non-BPS solitons in the same homotopy class. The 
energy of such objects is above the BPS bound, but they are static and topologically 
stable. 

To address this issue we study static solutions to the second order equations 
of motion. The kinetic energy is minimal for the shortest path in the complex 
plane, which means straight lines connecting equivalent points. In addition, the 
scalar potential has a saddle point at the origin in the target space and has ridges 
originating from this saddle point on the real axis. This means that there is a static 
real solution to the second order equation of motion for any value of L. For L = 1 
such a solution saturates the BPS bound. For L < 1, the kinetic energy dominates 
and the total energy is minimized by the straight line on the real axis in the complex 
<f) plane. Moreover, there are no BPS saturated solitons of the same homotopy class 
in this regime; the real soliton is stable because there is nothing to decay into. 
For L > 1, the kinetic energy does not dominate any more; there are other static 
solutions that are not straight lines in the complex plane that actually have lower 
total energy, the BPS saturated complex solitons. In this regime, the real soliton is 
unstable. We will discuss it in the next section. 

For L 1, the real non-BPS soliton is approximately given by 

4(x) = 2vr|, (24) 
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Figure 5: a) Energy of the real, non-BPS soliton as a function of the circumference 
L of the compact dimension; the BPS bound (dashed line), the analytical approxi- 
mation M so \ = U + K for small L (solid line), and the actual numerically calculated 
soliton mass (dots) are shown, b) Energy of the sphaleron, the real, unstable soliton 
for L > 1 as a function of the circumference L of the compact dimension; the BPS 
bound (dashed line), the linear analytical approximation for M sphal for large L (solid 
line), and the actual numerically calculated sphaleron mass (dots) are shown. 



where x ranges between —L/2 and L/2. For this configuration, the kinetic energy 
is given by 

A -./"V|) 2 ^, (25) 



L/2 



while the potential energy is 



U 



L ' 2 dx 167T 2 32L7T 2 

. L/2 X (2-cos27ra/L) 2 " 3^ ' 



(26) 



For L 1, the energy of the soliton therefore approaches M so \ = U + K. In Fig. |5]a) 
we show the BPS bound and the approximate soliton energy for small L, together 
with the numerically calculated energy of the non-BPS soliton. 



3.2 Sphaleron for L > 1. 

The real, static solution of the second order equation of motion persists even for 
L > 1. In this regime, the soliton is unstable and we will refer to it as a sphaleron, 
by analogy with the sphaleron in the Yang-Mills theory ||15|| . As in the Yang- 
Mills theory, the sphaleron mass in our problem will give the height of the barrier 
under which the solitons (/) tunnels into (r) and vice versa. Because the solution 
is real, the second order equation of motion with vanishing time derivative can be 
integrated. In fact, the equation of motion is identical to the equation describing the 
one-dimensional motion of a particle moving in the potential — V(x). The implicit 
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solution for d> is 



x-x = / ffl/y/V(6) + V , (27) 







where Vq is an integration constant (equivalent to the total energy of the particle). 
The solution <p(x) is periodic modulo 2ir with wavelength £(Vq). The constant Vq 
has to be adjusted so that the wavelength of the solution (which corresponds to the 
"time" it takes for the particle to move around the circle) is equal to the circumfer- 
ence L of the compact dimension, that is £{Vq) = L. This equation has one solution 
for any positive value of L. The method to determine Vq is similar to the procedure 
that was used to select I in the case of the BPS solitons. The energy of the real 
soliton/sphaleron is equal to 

M finha , = / d62JV(d>) + V n -V n L. (28) 



/7T 
d<j>2y/V(</>) + V Q -V L. 
-7T 



This energy can be explicitly determined in various limits. For Vq ^> 0, the real 
soliton of the previous section is obtained, with y/Vo = 2tx / L (this corresponds to a 
particle with so much kinetic energy that it hardly notices the potential), whereas 
for Vq = the solution is equal to the real BPS saturated soliton with L = 1 and 
energy equal to the BPS bound. When Vq is close to minus the minimum value of 
the potential, Vq ~ — 167r 2 /9, then L ^> 1 and the energy increases linearly in L, 

8tt 2 16tt. VS + 1 167r 2 r . . 

M sphal = 7f + — l Qg _ + __L + ..., (29) 

where the ellipsis indicate terms that vanish in the limit L — > oo. (This last situa- 
tion corresponds to a particle with just barely enough energy to reach the top of the 
hill). In Fig. ^b) we show the mass of the BPS saturated soliton and the approxi- 
mate sphaleron energy for large L in Eq. (p5|), together with the actual numerically 
calculated energy of the sphaleron. 



3.3 The Tunneling Action. 

For any value L > 1, there are four classical BPS saturated solitons, two in class Fx 
and one each in class T 2 and T3 . The two solitons in class Ti are marked by values 
of B that have the same magnitude but opposite sign. They are mapped onto each 
other in the complex <p plane by reflection in the real axis. In order to distinguish 
these two solitons, we will refer to them as the (I) soliton when B is negative, and 
the (r) soliton when B is positive. Tunneling mixes the two class Ti solitons, and, as 
a consequence, their mass is lifted above the BPS bound. The class T 2 and T3 BPS 
solitons do not mix; the tunneling action is infinite since the cycle would have to be 
moved across the poles (see Fig. f|); but the (I) and the (r) solitons can be deformed 
into each other without crossing a pole. The energy barrier that separates them is 
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therefore finite, see Fig. 5. In addition, at L ^> 1 the barrier is high, and the quasi- 
classical approximation is applicable. According to our previous considerations, the 
two shortened supermultiplets pair up to form a full representation, and the mass 
is lifted from the BPS bound. 



3.3.1 BPS bound on the tunneling action 

Here we will derive the BPS bound on the tunneling action for the specific model 
under consideration using the methodology outlined in Section |2|. The bound is 
saturated if the instanton configuration that interpolates between the (r) soliton at 
T = — oo and the (I) soliton at T = oo satisfies the two-dimensional BPS equation 
given in Eq. ([12]); in Section [3.3.2| we will use numerical methods to show that the 



instanton is indeed BPS saturated. 

In order to determine the surface terms in Eq. 
Euclidean action can be written as 



A 



BPS 



dt dx 



2V • S - (V A a) 



Tl), the BPS bound on the 



(30) 



where 



S 



Re(W) 
Im ( W) 



a = — 



lm(4>d x (, 
Im(0<9t<5 



(31) 



and we have used (V Aa) z as short-hand for d x a t — d t a x . Then application of Gauss' 
and Stoke's theorems converts the surface integral in Eq. ([50) into contour integrals 
over the boundaries of the surface, i.e. 



A 



BPS 



S ■ dn — d> a ■ dx 



(32) 



As noted in Section |2], this is the same equation derived in Ref. JT(J for the BPS 
bound on the energy of domain wall junctions. In the first integral in Eq. (j32j), 
dn is an infinitesimal vector perpendicular to the contour with length \dx\, pointing 
outwards from the enclosed area. In the second integral, dx is an infinitesimal vector 
tangential to the contour, and the contour must be followed counter-clockwise. 

The problem of calculating the Euclidean action is therefore equivalent to the 
calculation of the energy of a domain wall junctions, with the solitons corresponding 
to domain walls. The BPS bound on the action is completely specified by the 
boundary conditions. We have to deal properly with the fact that the x direction in 
our model is compact. In Fig. ^|we show the boundary conditions in the x,t plane; 
at t = —T/2 the field 4>(x,t) is equal to the (I) soliton, <fr(x, — T/2) = 4>i(x), where 
the (I) soliton is positioned so that it has its maximum energy density at x — 0. 
Similarly, at t = +T/2 the field 4>(x,t) is equal to the (r) soliton, <j>(x,+T/2) = 
4> r {x), where the (r) soliton is also positioned so that it reaches its maximum energy 
density at x = 0. We always have in mind the limit that T is very large. For 
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Figure 6: Boundary conditions on the world sheet. The initial (I) and final (r) 
solitons are located at the bottom and top sides of the rectangle at t = —T/2 and 
t = T/2, with the center of the energy density positioned at x = 0. Periodic (modulo 
2n) boundary conditions apply to the left and right side of the rectangle. The shaded 
area significantly contributes to the action. 



x = —L/2 and x = L/2, periodic (modulo 2tc) boundary conditions have to be 
imposed, </>(— L/2, i) = <ft(L/2,t) — 2n, as the space dimension is compact. 

The contour of the integrals in Eq. (|32|) follows the edges of the rectangle in 
Fig. |5|. Special care must be taken with the integrals over the vertical edges, at 
x = ±L/2. Naively, one might think that these contributions vanish, but because 
both the superpotential and the field are multi-valued, these integrals do in fact 
contribute. 

Before calculating the integrals in Eq. (|32|) , let us first determine the dominant 
contribution from the vertical wall in the large L limit. As stated in Section |^, this 
contribution is just L times the tension a of the horizontal "wall" in Fig. || If L 
is large, then the absolute value of B, the parameter that marks the (r) and (I) 
solitons, becomes large. In fact, the absolute value of B increases logarithmically 
with L. Except for points in space near the center of the soliton where the energy 
density is maximal, the soliton field is approximately equal to ±7r ± iB, where the 
second ± refers to the (r) and (I) soliton, respectively. For 5 = ±7r/2, there is a 
solution to the BPS equation that takes the form = ±tc + if(t), where f(t) — > ±00 
for t — > ±00. The tension of the horizontal "wall" a is equal to the tension of this 
solution. The leading contribution to AA, linear in L, is therefore 



AA = oh H 



32ttL 



V3 



arctanh 




) 



+ •••, 



(33) 
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where the ellipsis indicate subleading terms in L. We will now uncover these sublead- 
ing terms, which contribute to the pre-exponential factor in the tunneling amplitude. 
We will first determine the contributions from the spokes in Fig. |^ and then the hub. 



Spoke contributions 



We first simultaneously calculate the contributions of the first integral in Eq. (T32| 
over the left and right edges of the rectangle in Fig. |6], 



Ax 



T/2 

dt 

-T/2 



-L/2 



(34) 



In order to calculate Ai, it is necessary to exploit some symmetries of the soliton 
solutions. Apart from the translational invariance (modulo 2tt) 

<P(x + L) = 2tt + 0(x), (35) 

the soliton solutions have the following two symmetry properties 

Re [<f){L/2 + x)} = 2vr - Re [0(L/2 - x)} , (36) 

and 

Im [<j)(L/2 + x)] = Im [<j>(L/2 - x)] . (37) 

If the symmetry in Eq. (|36"D is preserved by the instanton configuration, the inter- 
polating field necessarily takes the form 



(f>(L/2,t) = ir + iB(L,t) 
<f>(-L/2,t) = -ir + iB(L,t), 



(38) 



at the left and right edges of the rectangle in Fig. ^, at x — ±L/2, where B(L, t) is a 
function interpolating between —B(L) at t — —T/2 and B(L) at t — T/2. However, 
from the identity 



W(±7T + IB) = 



8n 



V3 I 2 



7T 



± — h i arctanh 



1 i B 
—= tanh — 

V3 2 



(39) 



it is clear that the difference between the two integrands above does not depend on 
the function B(L,t) at all. The contribution A\ is therefore given by 



I671-2 



(40) 



This is just T times the BPS bound on the soliton mass, or the Euclidean action in 
the absence of the tunneling transition. 
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Next, we simultaneously calculate the contributions of the first integral in Eq. ( |32| 
over the top and bottom edges of the rectangle in Fig. 13, 



-L/2 

A 2 = 2 I dx 

-L/2 



Im(W)| t=T/2 -Im(>V)i t= 

-T/2 



(41) 



Here the integrands are just the constants of motion for the soliton solutions, so that 
A 2 = 2I r — 2Ii = 4I r , and I r and Ii mark the (r) and the (I) soliton, respectively. 
Therefore, in terms of B(L), the value of B that marks the (r) soliton for a compact 
dimension with circumference L, we obtain 



A, 



32tt r / 1 _ B(L) 
— — L arctanh —= tanh 



(42) 



For large L, such that B(L) ^> 2, this reduces to A 2 = aL, the dominant contribu- 
tion to the action that was already obtained in Eq. (|33|). 

Hub contributions 

We first simultaneously calculate the contribution of the second integral in Eq. 
over the vertical edges in Fig. || 



A, 



fT/2 






/ dt 


Im (c 


bdt$)\ 


/-T/2 







x=L/2 



Im 



\x=-L/2 



(43) 



Using the form of (p for x = ±L/2 given in Eq. (38), this contribution can be 
calculated and yields 



fT/2 

A 3 = -2n dtd t B(L,t) = -AttB(L). 

J -T/2 



(44) 



Finally, we simultaneously calculate the contribution of the second integral in 
Eq. ( |3~2"D over the horizontal edges in Fig. ||, 



L/2 



A4 = dx 

'-L/2 



Im (<f)d x (j)) I 



— Im (cj)d a 



(45) 



\t=-T/2 VT^XY) \t=T/2 

The (r) and (I) soliton are mapped onto each other by the transformation 

Im[(j)(x)} -> -lm[(f)(x)}. (46) 



At the top and bottom edges of the rectangle in Fig. || the field <ft therefore takes 
the form 

4>(x, T/2) = a(x) + ib(x), 
<fi(x, —T/2) = a(x) — ib(x), 



(47) 
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Figure 7: The (r) and (/) solitons and the instanton in the complex plane (compare 
Fig.(^) for the same configuration on the world sheet). The area S(L) is indicated 
by the shaded region. 



where a(x) and b(x) are real, and a(±L/2) = ±7r and 6(±L/2) = B(L). The integral 
can now be calculated and yields 

A A = -2 / dx {ad x b - bd x a) = AixB(L) - 2S(L), (48) 

J -L/2 

where 

/L/2 rx=L/2 pit 

dx bd x a = 2 dab = 2 dab, (49) 

-L/2 Jx=-L/2 J-ir 

is the area in the complex <p plane between the (/) and (r) soliton configurations 
(see Fig. 0). 

The hub contribution to the action is thus seen to be equal to minus twice the 
area traced by the instanton in the complex plane, Ah u b = A% + A A = —2S(L). In 
Ref. it was shown, within the context of domain wall junctions, that this result is 
valid in general and goes beyond the specific model we consider here. We also note 



that Ah u b is negative, in agreement with the general considerations in Ref. []T2| . The 
total semi-classical Euclidan tunneling action is 

AAbps = A 2 + A 3 + Ai, 

- 32n L arctanh ( -±= tanh - 2S(L). (50) 



We have calculated each contribution in terms of the functions B(L) and S(L). 
These functions depend only on the initial and the final soliton configuration. In 
the large L limit they take the form 

S(L) = 27rlogL + ---, (51) 
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and 

B(L) = logL + • • • , (52) 

where the ellipsis indicate terms that are finite or suppressed. We observe that the 
hub contribution depends logarithmically on L, in contrast with the situation for the 
domain wall junctions in models with infinite space dimensions, for example those 



considered in [LI], 0, where the hub contribution to the energy of the junction is 
finite. The difference is that for those models the energy density falls off exponen- 
tially fast away from the domain walls, whereas in the present model the energy 
density only falls off like the second power of the inverse distance to the wall. The 
circumference L of the compact dimension in a sense acts like an infrared regulator. 
In the large L limit, the dominant and subdominant terms in tunneling action are 
given by 

327T 1 

AA ~ —j= L arctanh— = — An log L + ■ ■ ■ , (53) 
v 3 v3 

where the subdominant term contributes to the pre-exponential factor in the tun- 
neling amplitude. 

3.3.2 BPS saturation of the instanton configuration 

The calculation in the previous section of the tunneling action hinges on the question 
whether a BPS saturated instanton configuration exists. The explicit form of such a 
configuration is not needed. In this section we use numerical analysis to address the 
question if a BPS saturated instanton configuration exists in this specific model. In 
order to numerically determine the instanton configuration, we embedded the model 
in d — 2 + 1 dimensions. One of the space dimensions is compact and the other 
space dimension represents the original Euclidean time. The new time is just an 
auxiliary construction. An initial configuration smoothly interpolates between the 
(I) and (r) solitons. The second order equations of motions are then evaluated and 
at the same time the system is cooled. The instanton configuration then emerges as 
a domain wall junction. 

We simulated the second order equations of motion on a lattice using a forward 



predicting algorithm similar to the one used in JIT], |12J . A complication arises due 
to the fact that both space and the target space have the topology of a cylinder, so 
that proper periodic boundary conditions had to be implemented. 

The lattice spacing was chosen to be much smaller than the width of the solitons. 
At the same time the size of the lattice was taken to be much larger than the soliton 
width. We performed the numerical calculation of the tunneling action for values 
of L ranging from L = 2 up to L = 6. The semi-classical tunneling action is 
only physically meaningful for larger values of L. However, the computation time 
increases rapidly with L. As shown in Fig. |8|, the calculated value for the tunneling 
action with the numerically determined instanton configuration saturates the BPS 
bound in Eq. (|50|). As there is no reason to think otherwise, we are confident that 
this saturation also occurs for larger values of L. 
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2 3 4 5 6 



Figure 8: Comparison between the BPS bound on the tunneling action from Eq. fl50p 
(solid line) and the actual semi-classical tunneling action from the numerical calcu- 
lation (dots). It is clear that the actual instanton configuration saturates the BPS 
bound. 



4 Quantum-mechanical description 

Under certain conditions on the superpotential W (and the Kahler potential /C, if it 
is non-trivial) it may be possible to develop a quantum-mechanical description of the 
tunneling of the distinct solitons resulting in the loss of their BPS saturation. Such 
an approach is valid provided there is a certain direction in the functional space 
which is much softer than all other "perpendicular" directions, and the tunneling 
occurs in this "soft" direction. Then all "perpendicular" degrees of freedom (there 
are infinitely many of them) can be integrated out and what is left is the quantum- 
mechanical motion of the center of mass plus the quantum mechanical dynamics of 
this particular soft degree of freedom. The motion of the center of mass is irrelevant, 
of course. It is always assumed that the soliton is at rest. 

The example considered in Section 3 contains no adjustable parameters (except 
L), so that one does not expect to find a specifically soft excitation mode in the 
soliton background. Therefore, in this case the quantum-mechanical description 
(with one bosonic variable) would describe the system only qualitatively. It is not 
difficult, however, to modify the superpotential to create a soft mode. To this end, 
let us consider the superpotential 

dW= n ^ d<S>, (54) 
2 - (3 cos $ 

which introduces a parameter j3 in the superpotential (0), without changing any- 
thing else. At very small j3 and large L, 

4 

log - > L > 1, (55) 

the x dependence of the imaginary part of becomes weak, and the softest mode 
corresponding to varying lm(0) is much softer than all other modes. (We treat L as 
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a fixed parameter, while (3 can be made as small as desired.) Then we can define our 
quantum mechanical variable B such that B is the value of lm(0) (for defmiteness, 
at Re(0) = ±7r). At the classical equilibrium 

B = ± \og(L/P) + const. (L > 1) . (56) 

The positive value of B corresponds to the (r) soliton, the negative to the (I) soliton. 
Both are classically BPS saturated. To roll over from positive to negative B, the 
system has to tunnel under the barrier at B = 0. 

Before passing to the discussion of supersymmetric quantum mechanics for the 
variable B we have to present a few formula from Section 3 in a modified form which 
includes the parameter 3, see Eq. (|54]). 

The superpotential with the additional parameter (3 takes the form 



W = — ; ^ arctan [ \ \ \ — ^ tan ^ 1 . (57) 



The poles of the scalar potential are located on the imaginary axis at 



so at small (3 they move further away from each other. The energy of the classical 
BPS saturated soliton is 

167T 2 

M BPS = Z = - -> 8vr 2 at (3 -> 0. (59) 

Finally, in the limit of large L the sphaleron energy is equal to 

M M 3271 ♦ V^HZL 167F 1 V /2 +^+V^ 

Jw sr> hai = J^bps ; arctan -= log . 

P V 7 ^ 7 ^ v 7 ^ 2 + /? s y2+?-^ 

16tt 2 

- 1^ 



(2 + /?) 2 " 

-> M BPS + 4tt 2 L + const, at /3 -»• 0. (60) 

We recall that the difference between the sphaleron energy and the classical soliton 
mass measures the height of the tunneling barrier. The tunneling action as a function 
of (3 becomes 



32vrL 2-3 A 
AA = arctanh. / £ -> 8%L log - at 8 -»• 0. (61) 

Now we pass to the quantum-mechanical treatment. To begin, let us recall some 
general aspects. A general algebraic consideration of the BPS soliton dynamics in 
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the extreme non-relativistic limit is carried out in [|nj. We will adapt it here for our 
purposes. 

The superalgebra we deal with has four supercharges and a central charge Z (the 
latter is real in the problem at hand). In the extreme nonrelativistic limit it can be 
represented as JTB] (in the soliton rest-frame) 



Ql = V2Zti®(X3 + ..., Q\ = V2Zt 2 <g> a 3 + ... , 

Q\ = I ® — L= [pa, + W\B)a 2 ] , Q\ = I <g> -±= \pa 2 - W'{B)a x ] , (62) 
l\ m 2\/m 



where 

p = —id/dB 

and m is an effective inertia coefficient for the variable B ("mass"). The dots in the 
first line in Eq. (p^) stand for higher order terms in the nonrelativistic expansion. 
Moreover, W (B) is the quantum-mechanical superpotential (not to be confused 
with the field-theoretic superpotential W). The algebra ( |6^) has, as its subalgebra, 
Witten's quantum mechanics [IT] . 



In principle, W(B) could be calculated from the underlying field theory, in the 



same way as it was done in a related problem in Ref. [16|. We will not do this 
calculation in full, since our purpose here is mainly illustrative. Instead, we will 
present a simplified superpotential which properly conveys qualitative features of 
the actual superpotential: (i) the existence of two zeros of W{B) at B = ±|£> | 
(double- well potential); (ii) the existence of the barrier at B near zero, with the 
maximum at zero, and with the height coinciding with M sp h a i — Z = Air 2 L. 
The corresponding superpotential and Hamiltonian have the form 

w = k(-^-+B^y 

H = f- + ^{W'f + ^-a 3 W" , (63) 
2m 2m 2m 

where 



B = log L/P 

log 4//3 
log 3 L/P 



k = forL , ° ' T , n (64) 



9 r log 2 4//? 

m = -L k 

2 log 2 L//3 

(Note that the overall two- by-two unit matrix following from the second line in 
Eq. (|6~2D is omitted in H. Its just replicates the states, in two copies each. We will 
keep it in mind.) 
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Bq is a function of L, which is logarithmically large and positive at large L. 
As L decreases, Bq decreases too, and goes to zero at L = 1. (To see that this is 
the case, one should inspect the exact expression for B rather than the asymptotic 
form in Eq . fl65|) valid at large L.) Below this point Bq < 0, and the Hamiltonian has 
no supersymmetric solution at the classical level. The Witten index of the system 
under consideration is zero. 

At positive Bq there are two classical solutions of the equation W = 0, cor- 
responding to two classical vacua. The tunneling between them was thoroughly 
studied, see e.g. [[0|. The one-instanton action is obviously 

4k 

Anst = [W(B = \B \) - W(B = -\B \)} = AW = — \B \ 3 . (65) 

The instanton transition is accompanied by a fermion zero mode, which suppresses 
all one-instanton amplitudes. The shift of the ground state from zero is given by the 
instanton-anti-instanton transition and is proportional to exp(— 2A- mst ) (for further 
details see ||18|| ). This can also clearly be seen from the Hamiltonian fl6~3f) which has 
a strictly conserved quantum number, [<j 3 , H] = 0. The one-instanton transition 
would flip the spin. The ground state of the Hamiltonian ( |63D is doubly degenerate, 
with one spin-up and one spin-down state of energy E ~ exp(— 2A inst ). If we include, 
in addition, the replication (which was mentioned above) due to the two-by-two unit 
matrix, we get that the overall number of states is four, as it must be in any non-BPS 
J\f = 2 supermultiplet. 

To conclude this section, we stress again that although the superpotential spec- 
ified in Eqs.(|63D and (|65|) correctly reproduces gross features of the tunneling phe- 
nomenon, it is definitely not the genuine superpotential that might arise in this 
problem should we decide to actually calculate it. In particular, it vanishes at small 
(3 only logarithmically, while actually one could expect a power-type suppression. 

5 Conclusions 

In a previous work, Ref. [[J, it was observed that in certain (1 + 1) -dimensional models 
with M = 2 supersymmetry and a compact space dimension, tunneling between two 
distinct, classically BPS saturated solitons lifts their mass above the BPS bound. 
The two shortened multiplets combine to form one regular M = 2 multiplet. One 
of the main new assertions in this paper is that the instanton that interpolates 
between the solitons in the Euclidean time is a 1/4 BPS saturated "domain wall 
junction." The semi-classical instanton action is determined by two central charges 
in the supersymmetry algebra. The action can therefore be determined solely from 
the solitons that are connected by the instanton; an explicit instanton solution is 
not necessary. 

Under certain conditions the description of the tunneling between the solitons 
in the field theory can be reduced to a non-relativistic supersymmetric quantum 
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mechanical model. This approach is valid if the direction in the functional space in 
which the tunneling takes place is much softer than all other directions (except for 
the zero modes), yet the tunneling action is large. The internal coordinate along 
this direction then forms the one remaining degree of freedom after all hard degrees 
of freedom have been "integrated out." 

We illustrated the above observations explicitly in a specific model. The semi- 
classical tunneling action was determined in the field theory. We used numerical 
methods to verify that the instanton indeed corresponds to the 1/4 BPS saturated 
domain wall junction. We also determined a range of parameters for which the 
tunneling can be described by a quantum mechanical model. Although possible in 
principle, we did not determine the actual superpotential of the quantum mechanical 
model in this range. Instead, we determined the parameters of a toy superpotential 
that is similar to the actual superpotential in all important features. 
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